Dynamical Diffraction Theory for Wave Packet Propagation in Deformed Crystals 
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We develop a theory for the trajectory of an x ray in the presence of a crystal deformation. A 
set of equations of motion for an x-ray wave packet including the dynamical diffraction is derived, 
taking into account the Berry phase as a correction to geometrical optics. The trajectory of the 
wave packet has a shift of the center position due to a crystal deformation. Remarkably, in the 
vicinity of the Bragg condition, the shift is enhanced by a factor (uj: frequency of an x ray, 
Aco: gap frequency induced by the Bragg reflection). Comparison with the conventional dynamical 
diffraction theory is also made. 

PACS numbers: 41.50.+h, 41.85.-p, 42.25.Bs 



X-ray diffraction has been used as a powerful tool 
for analyzing crystal structures. Conventional analy- 
sis of x-ray diffraction is based on a kinematical the- 
ory, namely the first-order Born approximation, and has 
played an important role in many fields beyond physics 
such as medical science. A theory incorporating the 
multiple-scattering process of x ray is called the dynami- 
cal diffraction theory Q , pioneered by Ewald 2] and Dar- 
win 0. In Darwin's theory for a perfect crystal, higher- 
order diffraction is taken into account by superposing 
the diffracted waves. A perfect crystal, however, suffers 
from an extinction effect; the intensity of the first-order 
diffracted wave is reduced by the higher-order waves. In 
order to avoid this difficulty, mosaicness is introduced 
into a crystal to suppress the contributions of the higher 
order diffracted waves by random interferences. Thus, a 
perfect crystal is not preferable for the structural anal- 
ysis, and applicability of the dynamical theory has been 
limited. However, recent technological advances for ex- 
tremely intense source of x ray such as a free-electron 
laser will open up new fields of research associated with 
the dynamical theory. 

On the other hand, for a deformed crystal, a dynamical 
theory was first developed dozens of years ago 
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Refs. 0,0, eigenequations for the Bragg reflected beams 
are translated into differential equations. In general, it 
is almost impossible to analytically solve such equations, 
and it is not so easy to extract a clear physical picture 
from the numerical solutions. We note a corresponding 
theory for electronic systems with dislocations pj |S| • 

These diffraction theories assume incident waves to be 
plane waves or spherical waves, and concern the propa- 
gation of their wave fronts. On the other hand, one can 
ask what happens if the incident wave is confined in the 
transverse direction, i.e., a narrow beam whose trajec- 
tory is well-defined. Naively one might expect that the 
trajectory of such a beam is always perpendicular to the 



wave front, which is expected in a wave propagation from 
Fer mat's principle. In this Letter, we reveal an anoma- 
lous feature of such electromagnetic beams beyond this 
naive expectation. We derive equations of motion for 
an x-ray wave packet in the presence of a crystal defor- 
mation. These equations of motion incorporate the dy- 
namical effect of a multiple scattering, and a geometrical 
phase associated with the wave dynamics, i.e., a Berry 
phase 9] , is taken into account in a natural way. We find 
that such a Berry phase in deformed crystals gives rise 
to a shift of the center position of an x-ray wave packet, 
and enhances the shift by several orders of magnitude in 
the vicinity of the Bragg condition. Such a gigantic en- 
hancement is physically related to a small group velocity 
near the Bragg condition. 

A wave packet has a finite width both in real and mo- 
mentum spaces due to the uncertainty principle. Such 
a wave packet feels "curvatures" in real and/or momen- 
tum spaces. A Berry phase describes these curvatures in 
real and momentum spaces, and plays a nontrivial role in 
the wave dynamics in constrained systems. For electro- 
magnetic waves, such a constraint is naturally realized; 
electromagnetic fields are always constrained to be per- 
pendicular to the propagation direction described by the 
wave vector k. The Berry phase for light appears for ex- 
ample in a twisted optical fiber in which the trajectory 
of the wave vector k makes a closed loop. In this case, 
the polarization plane rotates during propagation, and 
the rotation angle is represented by a Berry phase [To| . 
Another kind of constraint is realized in periodic sys- 
tems, where the dispersion relation has a band structure. 
When a physical state is constrained on each band, the 
Bloch function is associated with a Berry phase in mo- 
mentum (k) space. Namely, k space is curved for each 
band of the Bloch wave function, which influences the 
wave dynamics. Onoda et al. 11] constructed a general- 
ized geometrical optics including a Berry phase of light, 
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and derived a set of equations of motion for an optical 
wave packet. The Berry phase in k space gives rise to 
a transverse displacement of a wave packet, leading to 
the optical Hall effect. They were mainly concerned with 
photonic crystals for visible light, and investigated how 
to design the photonic band structure in k space. On the 
other hand, this Letter focuses on the x-ray diffraction; 
the main interest here is how to probe a spatial deforma- 
tion in real space. Therefore the relevant Berry phase is 
distinct from that discussed in 

Let us introduce a six-dimensional electromagnetic 
field vector; 
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where e and /z are the dielectric constant and the mag- 
netic permeability, respectively. Hereafter we put \x = 
h=l. The wave equations for F E ' H are 
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which are Hermitian. In these wave equations, we as- 
sume a periodic dielectric function in a crystal with re- 
spect to a primitive vector a: e(f + a) = e(f). In the 
dynamical theory, the eigenmodes of an electromagnetic 
wave are Bloch functions under the periodic condition. 
They are written in the six-dimensional representation 
as IV (r) = e lk ' r \Jt (r), where Ur (r) is a periodic 



function satisfying [Tg (r) = E/g (r + a), and z c rep- 
resents the polarization state. Figure da) shows the 
constant-frequency contour in k space, where the wave 
vector has a gap at the Bragg condition. Such a gap in 
k corresponds to an energy gap at the Bragg condition 
|fc| = \k— G|, as shown in Fig.^b). In the energy gap, the 
wave number becomes pure imaginary and no x ray can 
propagate through a crystal. Due to this gap, the disper- 
sion relation for the Bloch function has a band structure, 
and the Bloch functions are specified by a band index n; 
f n r (r). Using the Bloch functions, we construct a wave 
packet whose distribution function in k space has a sharp 
peak around the center of the momentum k c . The wave 
packet has a peak around the center position r = f c in 
real (r) space . We consider slowly varying perturbations 
which do not induce the interband transitions, and keep 
only the diagonal matrix elements with respect to the 
band index n. The distribution in k space gives an in- 
terference between components of different wave vectors, 
which can be expressed in terms of the Berry curvature 
in k space associated with each band. In addition, we in- 
troduce a crystal deformation whose spatial variation is 
small within the width of the wave packet. Such a slowly 
varying deformation leads to the Berry connection in r 
space, as we will show later. 
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FIG. 1: (a) Constant- frequency surface and corresponding 
wave vector and reciprocal vector k and (5, respectively, (b) 
Dispersion relation. At the Bragg condition, k and uj have 
gaps Ak and Auj. 



Let us derive equations of motion for such a wave 
packet. An effective Lagrangian for a wave packet is cal- 
culated in a similar manner as in |l2| to be 

L{r c -> k c , r c , k c , z c , z c ) = — uj n + k c • r c + ir c • (U Zc \^r c U Zc ) 
+ii-(U z jV£ c U Zc ) + i(U Zc \d t U Zc ). 

The first two terms on the right-hand side describe the 
wave front propagation, which exist even for plane waves 
or spherical waves. The last three terms represent the 
additional physics which appears only in the wave packet 
dynamics, leading to Berry curvatures. This effective 
Lagrangian is analogous to that for electronic systems 
|l2l H^. Variational derivative of the Lagrangian gives 
the following equations of motion; 
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where 



(fVfcW = tt ra k = {dr a U Zc \id k(3 Uz c ) ~ {dk p U Zc \idr a U Zc ), 

W AA , = ^Kxk c r^r c U nW%cf ), 
t^AA' = *(^nAfccfcl^fec nA'fc c r c )' 

\z c ) has two components representing the polarization 
states, and A is a polarization index. The vectors A-,k, 
called Berry connections in mathematics, are "vector po- 
tentials" in such a sense that their rotation gives the 
Berry curvatures. Each term in the above equations of 
motion can be interpreted as follows. The first terms on 
the right-hand side of Eqs. (1) and (2) reproduce geo- 
metrical optics. The second term in Eq. (1) represents 
effects of the crystal deformation, while the third one 
corresponds to a Lorentz force which is zero in a non- 
magnetic system. The second term in Eq. (2) represents 
the anomalous velocity of light 0? Q- giving rise to 
an optical Hall effect 11]. The third term represents 
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the displacement of a ray due to a crystal deformation. 
Equation (3) describes the time evolution of the polariza- 
tions. Since we consider a nonmagnetic system, £Vr = 0. 
We also assume that the perturbations are independent 
of time; Q t f — = 0. 

Consider a system with a deformation, as shown in 
Fig. [2j where the spatial variation of the deformation is 
small within the width of the wave packet. We assume 
that the wave packet has a sharp peak in k space around 
the center fc c , with its width much smaller than \G\. This 
implies that the peak in r space is much wider than the 
lattice constant. In this situation, the eigenfunction of an 
x ray is modified as T(r) — > T(f — u(r c )), where u(r) is 
a continuous function representing the atomic displace- 
ment 15]. In an equation of motion for fc Cj we assume 

= as the first approximation, resulting in k c = 0. 
Hence in the present level of approximation, does not 
affect the result. Thus we have only to consider the terms 
including and The Berry connection between 

the neighboring r c points is calculated as 



U z . 



: dU Zc 



d ~ 

U Zc (r-u(r c )) i — U Zc (r-u(r c )) 



where A T (k c ) is the Berry connection for a fixed polar- 
ization state in the absence of the deformation, and re- 
peated subscripts are summed over. Therefore, the Berry 
curvature term in the equation for r c is 



This is a product of the two factors. One is V? c u(r c ) due 
to the deformation, and the other is A r (k c ) originated 
from a band structure. Therefore, one can see that the 
Berry curvature represents an interplay between the 
curvatures in k and r spaces. Using these expressions, we 
can integrate the equation and obtain the displacement 
of an x ray due to the crystal deformation as 
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dm 
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is a group velocity, and fb = r c (t = 0). 



We show below that the factor enhances the dis- 

placement by several orders of magnitude in the vicinity 
of the Bragg condition. We assume that the wave vector 



I k z -(k-G) z 
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is in the vicinity of the Bragg condition, 
^j- ~ where Aco and Ak are gaps in co and fc, re- 
spectively. The Bloch function is expressed as f(r) c± 
Toe lk ' r + Tie l<K - k ~ G ^' r in two- wave approximation. From 




FIG. 2: Schematic illustration of the deformed crystal with an 
atomic displacement vector u(r), The trajectory of the center 
position of an x-ray wave packet is represented by the oriented 
curve. The amount of the shift can be as large as Ar c ~ 0.1 
mm. The broken line with an arrow represents the trajectory 
without the deformation. E< 
for s and p polarizations, respectively. 



and E p denote the electric fields 



this wave function, the enhancement factor - is cal- 
culated as 



dA 



r(3 
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where + and — denote the lower and the upper bands, 
respectively 0. Therefore, the solution to the equations 
of motion near the Bragg condition is 



vA±G 



G • 



(u(f c ) - u(fo)^j 



1 
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(4) 



This formula is the central result of this Letter. The di- 
rection for the shift is parallel to G, and the amount of 
the shift becomes maximum when G \\ u. Such a shift is 
enhanced by a factor -j^. This factor differs for polar- 
ization states, s- and p-polarizations, which are defined 
relative to the plane spanned by k and k — G. The dif- 
ferent size of the gap for each polarization state gives 
the different amount of the shift. When we inject an 
unpolarized beam, it splits into two beams with s- and 
p-polarizations, as shown in Fig. [21 

Now we give a realistic estimate for the enhancement 
factor and the shift Af c shown in Fig. [3 For an 
x ray, the dielectric constant can be approximated as 
e(r) = 1 — c^p/cj 2 , where cj p is a plasma frequency, and 
typically u^/ou 2 ~ 10 -6 [l(. Hence the value of the en- 
hancement factor for an x ray can be -j^ ~ 10 6 , which 
yields ^ ~ 10 -6 . For a deformation of \u\ ~ 0.1 nm, 
the displacement can be Af c ~ 0.1 mm. It means that 
the deformation in an atomic scale gives rise to a macro- 
scopic shift of the wave packet. Therefore, the factor 



4 



plays the role of a "lens" , which gigantically magni- 
fies the geometrical effect on the wave packet dynamics. 
Such a gigantic enhancement of the shift is a consequence 
of an energy gap and can be physically interpreted as fol- 
lows. In the vicinity of the Bragg condition, we have a 
very small group velocity v g . Such a slow wave packet 
has enough time to interact with electrons, resulting in a 
remarkably large optical effect. This kind of an enhance- 
ment at the edge of the energy gap is often studied in 
optical waves in photonic crystals Photonic crys- 

tals, however, have a large energy gap, and hence the 
enhancement is much smaller. 

We note that the distribution of the wave packet in k 
space is desired to be narrower than Ak. The components 
of wave numbers outside of this Ak width do not show 
the enhanced shift. To make the whole wave packet shift, 
an x ray with a wider spot size l r ~ (Ak) -1 ~ 0.1 mm is 
desirable, because its peak width in /c-space is as small 
as Ak ~ 10~ 6 \G\. The amount of the shift for such a 
beam is Ar c ~ l T . For a typical beam whose spot size is 
l v ~ 100 nm, only a small portion of the wave packet is 
shifted, and rest propagates without shifts. It is because 
the width in k space of this wave packet is larger than 
the gap Ak. It is then difficult to measure the shift. 

Let us compare the present theory with the conven- 
tional diffraction theory. In the conventional kinematical 
theory, only an equation of motion for k c is concerned. 
The equation for k c describes the propagation of the wave 
front. The term can deflect the beam; its amount, 
however, is not enhanced by the dynamical effect. The 
conventional dynamical theory deals with a part of the 
wave front which propagates through a deformation, and 
calculate the multiple scatterings which result in refrac- 
tion. Although a small part of the wave front has often 
been regarded as a wave packet, our theory reveals a 
novel physics of a wave packet, apart from wave front 
physics. The anomalous enhancement of the shift of the 
wave packet in this Letter occurs because the wave packet 
is not a monochromatic wave, but a mixture of waves 
with slightly different wave numbers. An interference 
between these waves gives rise to the enhancement dis- 
cussed in this Letter. This physics of a wave packet is 
absent in a plane wave or a spherical wave. 

It should be warned that such a displacement is com- 
pletely different from refraction. In refraction, the spot 
position on the screen is dependent on the distance be- 
tween a sample and a screen, because the refraction is 
caused by the shift of k c . On the other hand, the shift of 
r c obtained in this Letter is independent of the distance. 
Hence the shifts of k c and r c have different physical ori- 
gins and are experimentally distinguishable. 

In summary, we have developed a dynamical diffrac- 
tion theory for a wave packet propagation in a deformed 
crystal. Equations of motion for a wave packet are de- 
rived. The central result is the shift of the trajectory 
given in Eq. (4); the center of position of the wave packet 



is displaced by a crystal deformation parallel to the re- 
ciprocal lattice vector G, and the amount of the shift has 
maximum when G || u(f). Remarkably, such a shift of 
the wave packet is gigantically enhanced by the factor 
which is ~ 10 6 for an x ray. Since the factor is dif- 
ferent for each polarization state, the wave packets with 
no-polarization splits into two beams with independent 
polarizations. In contrast with the conventional dynam- 
ical theory which focuses on the wave front and requires 
numerically solving the Maxwell equations, our theory 
treats the narrow beam and provides an easy way to eval- 
uate a displacement of the wave packet for a given de- 
formation. An x-ray wave packet with a good coherency 
is useful to measure a deformation or dislocation in a 
crystal. 
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